Anisotropic non-perturbative zero modes for passively advected magnetic fields 
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A first analytic assessment of the role of anisotropic corrections to the isotropic anomalous scaling 
exponents is given for the d-dimensional kinematic magneto-hydrodynamics problem in the presence 
of a mean magnetic field. The velocity advecting the magnetic field changes very rapidly in time and 
scales with a positive exponent f . Inertial-range anisotropic contributions to the scaling exponents, 
Q, of second-order magnetic correlations are associated to zero modes and have been calculated 
non-perturbatively. For d = 3, the limit £ t-> yields ( 3 ■ = j - 2 + f; (2j 3 + f - 5j - 4)/[2(4j 2 - 1)] 
where j (j > 2) is the order in the Legendre polynomial decomposition applied to correlation 
functions. Conjectures on the fact that anisotropic components cannot change the isotropic threshold 
to the dynamo effect are also made. 



PACS number(s) : 47.27.Te, 47.27.-i 

Since Kolmogorov Q formulated his hypothesis, most 
of theories and models in turbulence have used as a 
key ingredient the restored local isotropy of small-scale 
structures, even in the presence of large-scale anisotropy. 
Dealing with such an idealized picture implies that the 
anisotropic effects, that almost every large-scale forcing 
indeed involve, are totally disregarded. 
Recently, some considerable effort has been done [|| [| to 
shed some light on the statistics of structure functions 
when taking anisotropy explicitly into account. When 
doing this, two major questions emerge: the first con- 
cerns the possibility of an universal nature of the scal- 
ing exponents of the separated isotropic and anisotropic 
contributions to structure functions; the second concerns 
the decay of anisotropic fluctuations, and consequently 
the validity of the local isotropy hypothesis. The state 
of the art on this subject, especially when looking at ex- 
perimental data, does not give unique answers. 
Here, our aim is to analyze, through non-perturbative 
calculations, the effects of anisotropy on anomalous 
(i.e. non-dimensional) scaling exponents of the mag- 
netic field correlations, within a kinematic magneto- 
hydrodynamics (MHD) problem (i.e. when reaction of the 
magnetic field on the velocity field is neglected). As the 
advecting velocity field that we consider is (^-correlated 
in time, analytical approach is possible: the main re- 
sult shown in this Letter is that the anomalous scaling 
exponent, Coj associated to the isotropic contribution is 
dominant with respect to the anisotropic ones, Q. In ad- 
dition, the entire set of anisotropic scaling exponents, Q, 
is given, showing the existence of a hierarchy related to 
the degree of anisotropy j, such that: Co < C2 < ■ • •• This 
result is the signature of the emergence of local isotropy 
at small scales. 

Such hierarchy relations hold analogous to those found 
in Ref. where the scalar advection in the presence of 
large-scale anisotropy is studied exploiting the field the- 
oretic renormalization group (RG). 

We also remark that the scenario here outlined is com- 



patible with those arising from the results shown in 
Refs. [|],[3| where, respectively, experimental and direct 
numerical simulation data of anisotropic turbulence have 
been analyzed. Considering the kinematic MHD prob- 
lem, the issue of the threshold for the appearance on the 
unbounded growth of the magnetic field is also present. 
We shall report a conjecture according to which we can 
argue that the anisotropic components do not play any 
effect in this sense. 

In the presence of a mean component B° (actually sup- 
posed varying on very large scales ~ L, the largest one 
in our problem) the kinematic MHD equations describing 
the evolution of the fluctuating part, B, of the magnetic 
field are §: 

d t B a +v-dB a = B ■ dv a + B° ■ dv a + Kd 2 B a (1) 
a = 1, • • • , d 

where the velocity, v, is a zero- mean Gaussian random 
process, homogeneous, isotropic and white-in-time, k is 
the magnetic diffusivity, and d is the dimension of the 
space. Both v and B are divergence free fields. The term 
B° ■ d v a in (Q) plays the same role as an external forcing 
driving the system and being also a source of anisotropy 
for the magnetic field statistics. 

The velocity is self-similar with the 2-point correlation 
function: 

(v a (r,t)v f3 (r',t'))=8(t-t') [d% - S a0 (r - r')] , (2) 

where S a p(r) is fixed by isotropy and scaling behavior, 
and scales with the exponent £, in the range < £ < 2: 



S a p(r) = Dr* (d + £-l)J a/3 -£ 



(3) 



The (^-correlation in time of v permits to exploit the 
Gaussian integration by parts (a comprehensive descrip- 
tion of this technique can be found, e.g., in Ref. ||) 
to obtain closed, exact equations for C a p(r,t) = 



1 



(B a (x,t)Bp(x + r,t)). After some manipulations of ([!]), 
such equations read: 



dt C a p — Sij didj C a p — (dj Sip) di C a j — 

(djS ai ) (d t C jP ) + (didjS^) {Cij + B°Bf) 
2nd 2 C a p. 



(4) 



A further equation for C a p follows from the solenoidal 
condition on B: 



d a C al3 = 0. 



(5) 



For what follows, it is worth emphasizing two properties 
of C a p: 

i) because of homogeneity, C a p is left invariant under the 
following set of transformations: 



and a < — > (3\ (6) 
as it follows from (||) after the 



ii) C a p{r) = C a p(-r) 
substitution r i— > — r. 
As shown in Ref. j|, in the isotropic case (i.e. B° = in 
our problem) anomalies appear already in the scaling ex- 
ponents of the second-order magnetic correlations, C a p, 
and have been calculated non-perturbatively by the au- 
thor. Anomalous scaling laws are associated with zero 
modes of the closed equations satisfied by the equal-time 
correlation functions. 

In Ref. , anomalous exponents for higher-order corre- 
lations have been calculated to the order £ by exploiting 
the RG technique. 

The extraction of anisotropic contributions to the 
isotropic scaling of C a p(r) found in Ref. ||, and the in- 
vestigation of their effect (if any) on the emergence of 
the dynamo effect are the main questions addressed in 
the present Letter. The main technical difference with 
respect to Ref. |J is that the angular structure of zero 
modes has now to be explicitly taken into account. 
In the presence of anisotropy, the most general expres- 
sion for the 2-point magnetic correlations, C a p(r), in 
the stationary state involves five (two in the isotropic 
case) functions depending on both r = \x — x'\ and 

* O " o 

z = cos 6 = B ■ r/r, where B is the unit vector cor- 
responding to the direction selected by the mean mag- 
netic field. Remark that the space is anisotropic but still 
homogeneous, so there is no explicit dependence on the 
points x, x' , but only on their distance. Namely, 

C a p{r) =J r 1 (r,z) S a p + T 2 {r 1 z)-^- + 

T 3 (r,z)^+T i (r,z)^^ + 
r r 

T 5 (r,z)B a B°p. (7) 

From the properties i) and ii) of C a p{r) one immediately 
obtains the following relations for the JF's: 



Ti(r,z) = J 7 i (r,-z) 

Fz(r,z) = ~^3(r, -z) 
Fz(r,z) = Ti{r,z). 



1,2,5 



(8) 
(9) 
(10) 



Substituting the expression (0) into (|) and using the 
chain rules, we obtain, after lengthy but straightforward 
algebra, the following four equations (corresponding to 
the projections over S a p, r a rp/r 2 7 B°rp/r and B° a Bp): 

[a ir 2 d 2 + hrd r + ci(l - z 2 )d 2 + d lZ d z + e^Ti + 
[fird r + gizd z + ji] T 2 + [fazrd r + hz 2 d z + m L z\T 3 + 
[o^p^Tr^iqi+nz^B 02 (11) 



[b 2 r 2 d 2 + c 2 rd r +d 2 (l- z 2 )d 2 + e 2 zd z + f 2 ] T 2 + 
92zT 3 + [k 2 + l 2 z 2 ] F h = (m 2 +n lZ 2 )B° 2 

a 3 d z Ti + b 3 d z T 2 + 

[c 3 r 2 d 2 + d 3 rd r + e 3 (l - z 2 )8 2 + f 3 zd z + g 3 ] T 3 + 
[j 3 z rd r + (k 3 + l 3 z 2 )d z + m 3 z] T 5 = n 3 B 



o 2 



a±d z T 3 + 

\b\r 2 d 2 + c^rdr + d 4 (l — z 2 )d 2 + e 4 z<9 z + / 4 ] T§ 
94B° 2 , 



(12) 



(13) 



(14) 



where the coefficients dj, 6j, • • • fj are cumbersome func- 
tions of £ and d and will not be here reported for the 
sake of brevity. Without loss of generality, we have fixed 
D = 1 in (|J), and we have neglected all terms involving 
the magnetic diffusivity k, our attention being indeed fo- 
cused in the inertial range of scales, i.e. r\ <C r <C L where 
rj = k 1 /^ is the dissipative scale for the problem. 
With the substitution of the expression (g), the solenoidal 
condition (|J) splits into the following couple of equations: 



[r9 r + (d-l)]^i 
\zrd r + d 7 — z 2 d z 



[rd r - 
A F3 



zd z ] T 2 
= 



(15) 



d z T 2 + [rd r +d]T 3 + [z rd r + (1 - z 2 )8 z ] T 5 = (16) 

associated to the projections over rp/r and respec- 
tively. 

From the relation (|h]) and Eqs. (|l^) and (|l6|) it then fol- 
lows that only two functions, J-'s, in (Q) arc independent. 
According to the old idea of Kolmogorov, in cascade-like 
mechanisms of transfer of energy towards small scales, 
anisotropy present at the integral scale should eventually 
decay during the (chaotic) transfer. One could thus argue 
that, at least at small scales, anisotropic corrections to 
the isotropic contribution would be smaller and smaller 
as the order of anisotropic contributions increases. For 
Navier-Stokcs turbulence in channel flow, such a picture 
has recently been confirmed by Arad et al in Ref. || . 
As we shall see, the above physical hint can be easily ex- 
ploited if one decomposes functions J-'s on the Legendre 
polynomial basis. Accordingly, we have: 



^(r,z) = £/f (r)P 2l (z) 
1=0 



1,2,5 



(17) 



2 



(18) 



;=o 



where the separation of even and odd orders in ( |17| ) and 
( |l8| ) arises as a consequence of the symmetries expressed 
by the relations (|J) and (Q), respectively. As larger Vs 
correspond to higher order anisotropic contributions, we 
thus expect that, when scaling behavior sets in (i.e. for 
f] <C r <C L), we shall have: 



Ji (r) cx r"i 



with 



(19) 



We would like to obtain equations for the /, (r), to be 
then solved for the Q 1 ': to do this, in the Eqs. (O)-(ll^) 
we have to express quantities like z^d™ (j,m — 0,1,2, 
with j 1 and m =/= 2) in terms of Legendre polynomials. 
Recalling to this purpose well-known relations involving 
these latter (see, e.g., Ref. Q), we obtain general ex- 
pressions as e.g.: 



1=0 



(2l + l)J2f£ + i+i(r) 



9=0 



(20) 



from which we notice that an arbitrary Z-order 
anisotropic contribution is coupled to all larger orders. 
The resulting equations arising from are thus 

not closed. 

Closed equations for the f's can actually be obtained by 
exploiting (|l9|), i.e. by using the hypothesis of a hierar- 
chy in the self-similar behavior of the f's. Accordingly, 
in Eqs. ([ll])-([l6]) at a given order j, for each function 

fl we need to retain only its lower order contributions 
with I < j. It is worth noticing that we can control the 
validity of this (physical) assumption in a self-consistent 
way, at the end of our calculation. 

As a result, one obtains the following (infinite) set of 
closed differential equations, valid for j even: 



p'(2) 



f(2) 



air 2 /-' (1) + hrf'/V + Cl /j 1} + dirff + erf?' + 
girf^l + hff\ + kiff} 2 = B° 2 h 6 j2 

f (l) . , 2t" ( 2 ) i /( 2 ) i J A 2 ) i 

aafj + hr fj + c 2 rf> + d 2 f) + e 2 f)A l + 
g 2 ff\ = B o2 ]2 5 : 



J2 



,(1) . , A 2 ) , 2,"(3) . . /(3) , j;(3) 

a 3 fj + hf) + c 3 r f^-> + d 3 rf^_{ + e 3 />j' 1 + 
9 3 rf'£l+h4%=B° 2 k 1 S j2 



(21) 



(22) 



(23) 



aiff} x + hr 2 /^ + Ci rf'£l + d t ff% = B° 2 e 4 5 j2 (24) 



„2/'(5) 



f'(5) 



f(5) 



2j 



3 — rf' (3) - 



2J-1 



(25) 



(2j 

(i-2)(j 



d ff-i 



2j-3 



/(5) 



2j 



A5) 
1 3 -2 



0, 



(26) 



where coefficients <Zj, 6j, • • • (different from those defined 
in (pd])-(p"4|)) depend only on £ and d. For j = 0, the co- 
efficients relative to and /._ 2 (and their derivatives) 

are zero and the resulting equations for f^ and are 
exactly as in Ref. |J. The structure of the above equa- 
tions fixes the relation between the scaling exponents rel- 
ative to different /'s. Indeed, when searching for power 

law solutions /j'' (r) cx r**i , in order to obtain balanced 
equations the 'oblique' relations must hold: 



(i) 



(2) 



A3) 



C 



(5) 
i-2- 



(27) 



We are now ready to show that nontrivial scaling behav- 
iors for the /'s take place due to zero modes, i.e. the 
solutions of the homogeneous problem associated to 
Eqs. (pl])-(p6|). To that purpose, let us consider such dif- 
ferential problem with no forcing (i.e. B° = 0). If, when 
looking for power law solutions, we exploit ( ]27| ) and the 
fact that only two functions of the f's are independent, 
our differential problem is mapped into an algebraic one. 
In this latter the emergence of zero modes reduces to 
imposing the existence of non-zero solutions of a 2 x 2 
homogeneous linear system. The calculation is lengthy 
but straightforward and the following expressions for the 
zero-mode exponents are obtained 113] : 



C 



1 



± _ 

2(d 

6d 3 + 4£ 2 d 



— {2£ + d 2 -d- [-2d 3 £-2d 2 £ 2 - 



+ 10 d £ + 20 d j - 20 d - 8 £ - 
8j + 4d 2 j 2 + 2£ 2 - 4£ j 2 + 17 d 2 -8dj 2 + 
8^j + Ad 3 j + Ad 2 j^ + Adj 2 ^ + Aj 2 -l6d 2 j - 

12d£j + d 4 ± 2VK(d- 1) (2-f) } (28) 



where 

K = (d- 1) {d 3 + Ad 2 j- 5d 2 + 2d 2 £ + i 2 d + Ad(,j - 
6d£ + 8d- I2dj + Adj 2 -£ 2 +4^ + 8j-80' - 
4-4j 2 + 4£j 2 )- 

Some remarks are in order. First, Q$ coincides with the 
isotropic solution obtained by Vergassola in Ref. 0] , the 
admissibility of which has been proved by the author. 
Second, £q~ diverges as r~ 2 at the dissipative scale 77. 
Exponent is thus not admissible. Third, > QJ 
for all j. This means that, in the inertial range of scales 
(i.e. r/L <C 1) the leading zero-mode solutions for j > 2 
are associated to C~ . We can thus define the leading set 

of zero modes Q as: Co = Co" i = Cj 7 ^ or — 2- 



3 



In particular, for j = 2, the asymptotic limits £ -C 1 and 
d^> 1 are, respectively: 



<2 



2£ 



(d-l)(d + 2) 



+ 0(£ 2 ); C2 



2£ 
d 2 



0(l/d 3 ). (29) 



Let us briefly discuss the infrared (IR) behavior of Q . In 
the absence of forcing terms, there is no way to satisfy 
the IR boundary condition (i.e. C a p(r) = for r i— ► oo): 
Q indeed diverges for r h oo. As a consequence, zero 
modes for j > 2 are not globally acceptable. For j = 2 
the situation changes completely: in this case Eqs. (^l|)- 
( p4| ) are forced and, as in Ref. IR boundary con- 
ditions can be satisfied by matching at the large scale 
L zero-mode solutions with those of the inhomogeneous 
problem. From the above considerations, it also follows 
that zero modes associated to Q become acceptable for 
all orders j when a fully anisotropic forcing term (i.e. pro- 
jecting on all Legendre polynomials) is added in the right 
hand side of (|l|). 

Finally, autoconsistency of our solution for Q, that is 
the validity of the hierarchy in ( |T^ ) , can be immediately 
checked from Fig. 1 where the behaviors of few Q as 
functions of £ are shown in the three dimensional case 
for j = 0,2,4 and 6 (from below to above). It is easy 
to verify that the increasing of scaling exponents with j 
actually holds for all values of j and d. 



condition as a criterion to select the emergence of an 
unbounded growth, we can conclude that there is no ef- 
fect played by the anisotropic components. Indeed, it is 
easily verified from (]2q ) that, for all d's, Q is real for 
£ G [0,1]. In conclusion, we have presented a system 

where the extraction of anisotropic contributions to the 
anomalous scaling of the equal-time magnetic correlation 
functions can be performed in a nonperturbative way. We 
have calculated the entire set of universal anomalous ex- 
ponents, Cji an d we have given an analytic assessment 
of the dominance of the fundamental exponent associ- 
ated to the isotropic shell. More generally, the hierarchy 
Co < Qi < ■ ■ ■ Cj ' < ' • • has been proved. The picture here 
drawn is in agreement with recent findings by Antonov in 
Ref. B, where the passive scalar problem is studied, and 
by Arad et al in Ref. M for Navier-Stokes turbulence. 
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FIG. 1. Behavior of Q vs £ for d = 3 and (from below to 
above) j — (heavy line) j = 2,4 and 6 (thin lines). 

The expression for Q allows us to make some con- 
jectures on the role played by anisotropic effects on the 
emergence of the dynamo effect. It is known ||,|l3| that in 
the isotropic case an unbounded growth of the magnetic 
field (dynamo effect) arises for £ > 1. The question ad- 
dressed here is whether anisotropic component can con- 
tribute to destabilize the system, shifting toward smaller 
values of £ the threshold for the dynamo. We note that, 
in the isotropic case, dynamo arises when the exponent 
related to the admissible zero mode becomes complex. 
In this case, zero-mode solutions have sinusoidal com- 
ponents, a fact that makes possible their matching with 
the appropriate boundary conditions also in the absence 
of forcing (i.e. the system is self-maintained). This hap- 
pens for £ > 1, £ = 1 being the threshold. Taking such 
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